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4.7 Dispersion in an oscillatory shear flow 

Relevant to the convective diffusion of salt and/or pollutants in a tidal channel, and chemicals 
in a blood vessel, Let us examine the Taylor dispersion in an oscillating flow in a pipe. Let 
the velocity profile be given, 

u = U s (r) + 3? [U w (r)e- iujt ] , < r < a. (4.7.1) 

The transport equation for the concentration of a solvent is recalled 

dC_ 8{uC) _ (d 2 C }_d_( r dC\\ (472) 
dt dx \ dx 2 r dr \ dr J J 

Assume the pipe to be so small that diffusion affects the whole radius within one period or 
so, i.e., 

2?r a 2 ,„„„x 

To 7T (4.7.3) 

to U 

We shall be interested in longitudinal diffusion across L much greater than a. Let U Q be the 
scale of U and 

x = Lx', r = ar', u = U u , t = ^t', Q = (4.7.4) 
Equation (4.7.2) is nomalized to 

dC Ua a dju'C) _ a 2 d 2 C I d_ ( dC\ 
~W + ~DL dx 1 ~ U ~dx^ + ~r~dr \~dr ) 

Let the Peclet number Pe = Ua/D = 0(a/L) be of (4.7.5) becomes 

dC | cPc d(u'C) =c2 d 2 C | 1 8 ( r ,dC'\ 
dt' dx' dx' 2 r dr' \ dr' J 

with the boundary conditons 

dC 

— = 0, r' = 0,l (4.7.7) 

with 

u' = U' S + m' w e~ lM (4.7.8) 
For brevity we drop the primes from now on. 
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1 ( 4 - 7 - 9 ) 



4.7.1 Multiple scale analysis-homogenization 

For convenience let us repeat the perturbation arguments of the last section. 

There are three time scales : diffusion time across a, convection time across L, and 
diffusion time across L. Their ratios are : 

t. L i I 1 

D 'U~ '~D~ '~e'7 2 

the smallest time scale being comparable to the oscillation period. Upon introducing the 
multiple time coordinates 

t, ti = et, t 2 = e 2 t (4.7. fO) 

and the multiple scale expansions. 

C = C + eC 1 + e 2 C 2 + .... (4.7.11) 

where Cj = Ci(x, r, t, t±, t 2 ), then the perturbation problems are 
0(e°): 

dt r dr ( <9r ) ( 7 ) 

with the boundary conditions: 

^ = 0, r = 0,l. (4.7.13) 

dti dt dx r dr \ dr J 

dC x 



0(e): 



with: 



dr 



= 0, r = 0,l. (4.7.15) 



0(e 2 ): 



with 



^£ + ?£l + ^ + p e ^ll = ^ + I A (4.7.16) 

<9t 2 9ti dt dx dx 2 r dr \ dr J 

dC 

^ = 0, r = 0,l. (4.7.17) 

Ignoring the transient that dies out quickly and focusing attention to the long-time 
evolution, i.e., t\ = 0(1), the solution at O(e ) is 1 

C = C (a;,t 1 ,t 2 ), (4.7.18) 



1 Strictly speaking the solution is 

oo 

Co = C 00 {x, ti,* 2 ) + ^ C „(x, t u t 2 )e- { < )2t J (k' n r) 
o 

where k' n is the n— th root of Jg(fca) = 0. The series terms die out quicly in t 3> 1 and ti <C 1 , leaving the 
limit Coo which is independent of t. (Dr. E. Qian,1993) 
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At 0(e), let the known velocity be 

u = U s {y) + ^(u w {y)e- int ) 

then 

dC dd 



+ ^ + Pe{U s + ^[U(r)e^ t }} 



dU dt 

Denoting the period average by overbars, 



dC _ 1 d ( dCi 
dx r dr \ dr 



f = rJ, ' 



dt 



and taking the period average, 

dC 



+ PeU. 



dC Id ( dd' 



with 



dti dx rdr\ dr 

dd 



dr 



= 0, r = 0, 1 



Let us now integrate (or average ) across the pipe, and get 

where angle brackets denote averaging over the cross section. 

1 f 1 



(h) 



7T JO 



2nrh dr 



Now subtract (4.7.23) from (4.7.20) 
dd 



dt 



+ Pe {U s + 9fc [U w e tnt ] } 



where 



dC _ 1 d ( dd 
dx r dr \ dr 

U = U a (y)-{U 9 ) 



(4.7.19) 
(4.7.20) 



(4.7.21) 



(4.7.22) 



(4.7.23) 



(4.7.24) 



(4.7.25) 



is the velocity nonuniformity 

Now d is governed by a linear equation, we can assume the solution to be proportional 
to the forcing and composed of a steady part and a time harmonic part, i.e., 



then 



d = Pe^ {B s (r) + 3? [B w (r)e- iQt ] } 
1 d ( dB? 



r dr \ dr 



U(r) 



(4.7.26) 
(4.7.27) 
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and 



r dr \ dr j 
with the boundary conditions 

dB, 

and 



1 d l r ^L\ +Ww = U w {r) 



dB n 



dr dr 
After solving for B S ,B W we go to 0(e 2 ), i.e., (4.7.16) 



0, r = 0,1. 



dC dC 1 + dC l 
dt 2 dt 1 dt 

+ Pe 2 {(U s ) + U s + ^ 

d 2 C 1 d ( dC 2 
+ 



U w e 



-int 



}{B s + U[B w (y)e- int ]} 



d 2 C 
dx 2 



dx 2 r dr y dr J 
which is a linear PDE for C 2 . From(4.7.26) and (4.7.23) we find 



= ~Pe 2 ^(U a ) {B s (r) + M[B w (r)e-™]} 



dx 2 



It follows that 



dt 2 dt 

+ Pe 2 {U s + ^[U w e^ int ]}{B s + ^ 

d 2 C | 1 d fdC 2 \ 
dx 2 r dr y <9r y 

Taking the time average over a period, 



B w (r)e- int ]} 



d 2 C 
dx 2 



dC 
dh 



<9 2 C <9 2 Q) 1 9 / dC 2 
+ -— I r- 



with 



' dx 2 dx 2 r dr y <9r 



<9r 



= r = 0,l 



Averaging (4.7.33) across the pipe, we get 



dC E d 2 C 



du 



dx 2 



with 



E = 1 



Pe 2 {(U s B s ) + ^(U w B* w )} 



(4.7.28) 



(4.7.29) 



(4.7.30) 



(4.7.31) 



(4.7.32) 

(4.7.33) 
(4.7.34) 

(4.7.35) 
(4.7.36) 
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which is the effective diffusion coefficient or the dispersion coefficient. The first part is of 
molecular origin; the second part is due to fluid shear. 
Finally we add (4.7.23) and (4.7.35) to get: 



This describes the convective diffusion of the area averaged concentration, which is cer- 
tainly of practical value. 

After the perturbation analysis is complete, there is no need to use multiple scales; we 
may now write 



still in dimensionless form. This equation governs the convective diffusion of the cross- 
sectional average, after the initial transient is smoothed out. 

Homework: Find the dispersion coefficient E in the oscillatory flow in a circular pipe 
and carry out the necesary numerical calculations. 

Homework (mini research) : Find the dispersion coefficient E in the oscillatory flow in 
a blood vessel with elastic wall. 




(4.7.37) 




(4.7.38) 



